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Recent experiments report that the long looked for thermotropic biaxial nematic phase has been 
finally detected in some thermotropic liquid crystalline systems. Inspired by these experimental 
observations we concentrate on some elementary theoretical issues concerned with the classical 
sixth-order Landau-deGennes free energy expansion in terms of the symmetric and traceless tensor 
order parameter Q a f3- In particular, we fully explore the stability of the biaxial nematic phase giving 
analytical solutions for all distinct classes of the phase diagrams that theory allows. This includes 
diagrams with triple- and (tri-)critical points and with multiple (reentrant) biaxial- and uniaxial 
phase transitions. A brief comparison with predictions of existing molecular theories is also given. 
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I. INTRODUCTION 

The biaxial nematic phase, predicted theoretically by 
Freiser 0, Q over 35 years ago, is one of the peren- 
nially challenging problems of experimental soft-matter 
physics. Although discovery of this phase was made by 
Saupe and co-workers in a fine-tuned lyotropic liquid 
crystal system in 1980 ;:| only in the past three years- 
and following several earlier attempts that proved un- 
successful in this regard (for a comprehensive review see 
e.g. 0, Q)- strong experimental evidence has become 
available that this phase can also be made stable in ther- 
motropic liquid crystalline materials @, 0, B This 
discovery raises the emerging theoretical problem of what 
mechanism is responsible for the stability of thermotropic 
biaxial nematics, especially for bent-core systems @, 0] 
and for tetrapode-like molecules [1, Qlj , where this phase 
was shown to be stable. 

There are two nematic phases of distinct symmetries. 
The ubiquitous uniaxial nematic phase has the 2?oo/i 
point group symmetry [ll], [l2|, [H, [l4| , which results in 
the definition of a single mesoscopic direction, known as 
the director. The director is a unit vector, denoted n, 
with the directions n and — n being equivalent. One con- 
sequence of this is that there are two different principal 
components of a second rank tensorial property, such as 
e.g the magnetic susceptibility. Generally, two uniax- 
ial nematic phases are distinguished: prolate (Nu + ) and 
oblate (Njj_). The prolate uniaxial states usually occur 
for rod-like molecules while disc-like molecules yield the 
oblate uniaxial states. As opposed to the uniaxial ne- 
matic phase, the biaxial nematic phase, denoted Nb, is 
characterized by three orthonormal directors, the Gold- 



stone modes, which we denote {1, m, n = 1 x rh}. Due 
to overall lack of polarity of the known biaxial nemat- 
ics, one finds that 1 and —1, rh and — m and n and — n 
directions are equivalent. That is, from the symmetry 
point of view the biaxial nematic phase is a structure of 
T^2h point-group symmetry and the corresponding sec- 
ond rank tensorial property has three different principal 
components. 

Generally, first- and second order phase transitions are 
observed experimentally between the isotropic phase and 
different nematic phases and between the nematic phases. 
The phase sequence of I so «-> (Ny-) «-> (Nb) <-> (Njj+) 
<-> (Nb) <-> (Nu-) <-> (I so) is found with decreasing tem- 
perature where the brackets indicate that 
some of the phases may not appear. In particular, the 
amazing reentrant uniaxial and isotropic phases are ob- 
served in lyotropic systems (see e.g. jl6| | and references 
therein) . 

On the theoretical level, possible effects of molecular 
structure on nematic order have been studied. More 
specifically, molecular field theories of single-component 
systems consisting of biaxial molecules and interacting 
via hard-core or continuous potentials were shown to p ro- 
duce a stable biaxial phase [J OS H H HI HI HI] • A 
similar scenario emerges from computer simulation stud- 
ies [2^ [2jj25U2^, H3, SI, HI] and from Landau treatments 
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Out of the theories cited the simplest description of 
the uniaxial and biaxial nematic phases is one offered 
by a sixth-order Landau-deGennes free energy expansion 
in terms of the alignment tensor Q a p ©. The theory is 
generally employed to interpret experimental data as well 
as to classify possible topologies of the phase diagrams. 
Therefore it seems quite important to know, if possible, 
an analytical form of all distinct classes of the phase dia- 
grams and limitations on them that can be derived from 
this simple theory. This task has on ly p artly been re- 
alized so far 0, [lj, HI HI, M, H, m- Nonc of the 
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papers cited shows, however, a full spectrum of predic- 
tions of this theory. The closest to the ideal is the paper 
by Prostakov but also there not all cases/analytical 
solutions have been given. 

Owing to the current excitement in the field of ther- 
motropic biaxial nematics we think it is important to 
re-examine this fundamental theory. We give analytical 
formulas for all distinct classes of the phase diagrams 
the model can predict and for their stability range. We 
hope this will be of some help for experimentalists in an- 
alyzing experimental data on biaxial nematics and will 
bring partial order to existing molecular predictions on 
this phase. 

This paper is organized as follows. After a brief discus- 
sion of the Landau-deGennes theory in Sec. II, we give 
analytical solutions for the phase diagrams in Sec. III. 
Section IV is devoted to a short discussion. 
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FIG. 1: Allowed variation of the independent degrees of free- 
dom Tr(Q 2 ) and Tr(Q 3 ) (shaded area) and identification of 
the corresponding phases. Shown are also lines of constant 
biaxiality parameter, uj. 



II. LANDAU-DEGENNES FREE ENERGY 

The best way to account for a symmetry change that 
takes place across a phase transition is by referring to an 
order parameter. For a phenomenological description of 
the nematics the relevant order parameters are tensors 
built out of the directors. Among these the leading order 
parameter is the second rank symmetric and traceless 
alignment tensor Q. In a standard parametrization Q 
can be written as 

Q = -^L (3n<g>n- 1) + m® m) , (1) 

where the directors {1, m, n} are identified with eigen- 
vectors of Q corresponding to the eigenvalues Ai = 

-75 + 75 ' A2 - "75 " % A3 = - Al - A2 = V§ 90 ' rc " 
spectively. The parametrization |T|) for Q is chosen such 

that the formula for F is kept concise. The isotropic state 
is stabilized when all three eigenvalues of Q arc equal and 
hence vanish, which yields Q = 0. For the T>ooh- sym- 
metric uniaxial states two out of the three eigenvalues of 
Q are equal, i.e., q Q ^ 0, q 2 = or q ^ 0, q 2 = V3qo 
or q 7^ 0, q 2 = — v3(?o- In the general case, Q has three 
different real eigenvalues that account for the P2/1- sym- 
metric biaxial state. A microscopic interpretation of the 
alig nment tensor for simple molecular models is found in 
[l8t [35| and can easily be extended to the more general 
cases. 

The Landau-deGennes phenomenological theory of 
non-chiral systems is implicitly based on the hypoth- 
esis that equilibrium properties of the system can be 
found from a non-equilibrium free energy, constructed 
as an C(3)-symmetric expansion in powers of Q. The 
only restriction on the expansion is that it must be sta- 
ble against an unlimited growth of the order parameter. 
There are two types of 0(3) invariants that can be con- 
structed out of Q, which involve traces and determinants 
of powers of Q. But determinants can be expressed in 



terms of traces and all traces of Q" with n > 4 are poly- 
nomials of Tr(Q 2 ) and Tr(Q 3 ) In addition, Tr(Q ) 
and Tr(Q 3 ) are bounded by the inequality 

^Tr(Q 2 ) 3 - Tr(Q 3 ) 2 = 

i (Ax - A 2 ) 2 (2Ai + A 2 ) 2 (Ai + 2A 2 ) 2 > 0, (2) 

which is fulfilled as equality for the uniaxial phases. 

A coordinate-independent form of the inequality (|2|) 
is obtained by a very convenient re-parametrization in 
Tr(Q 2 ) and Tr(Q 3 ) that uses just two scalar parameters: 
q and < u) < 1. They are introduced through the 
relations 

Tr(Q 2 ) = q 2 = q 2 a + q\ = I 2 (3) 

y/6 Tr(Q 3 ) = q 3 (l - u) = q% - 3q Q ql = I 3 , (4) 

where \q\ is the norm of Q and u> serves as a normalized 
measure of phase biaxiality. The T> 2 h- symmetric biaxial 
state is characterized by lu > with maximal biaxiality 
being accomplished for w = 1. For the uniaxial phases 
w = 0. In addition, for uniaxial Q-tcnsors a transforma- 
tion u — > [{Qafj — cS a f3)u a U/3], where c is an arbitrary 
constant making the bilinear form [...] positive-definite, 
transforms a unit sphere |u| = 1 into an axially sym- 
metric, prolate- (q > 0) or oblate (q < 0) closed surface. 
Hence the sign of <j, being consistent with the sign of 
Tr(Q 3 ), allows one to distinguish between Nu + (q > 0) 
and Nu- (q < 0) phases. Actually q and to can serve as 
invariant measures of order in uniaxial (q 0, lu = 0)- 
and biaxial (q ^ 0, to 7^ 0) nematics. For the isotropic 
phase g = 0. The allowed variation of Tr(Q 2 ) and Tr(Q 3 ) 
and consequently also of q and u>, along with the identi- 
fication of different nematic phases, is shown in Fig. |T]). 

In the absence of electric and magnetic fields the bulk 
free energy for the isotropic- and the nematic phases has 
the form 

F(Q) - F[Tr(Q 2 ),Tr(Q 3 )] = F[I 2: I 3 ] = F[q,u]- (5) 
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The minimal coupling Landau expansion of F that ac- 
counts for the biaxial nematic phase has to be taken up to 
6th order with respect to Q. This theory, also known as 
Landau-deGennes free energy of biaxial nematics, reads 
(see e.g. 
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(6) 

(7) 

(8) 



In Eq. ([6]) the F Q -part represents the unimportant free 
energy of the reference isotropic phase; F u is the free en- 
ergy of the uniaxial phases (cj = 0) and the remaining 
two terms represent biaxial contributions. The coeffi- 
cients of the expansion generally depend on temperature 
(inverse density) and other thermodynamic (control) pa- 
rameters. In what follows we will only keep an explicit 
dependence on the temperature. In particular, the coef- 
ficient a = a (T — T*) with T being the absolute temper- 
ature, is the only term in the expansion that is assumed 
to be explicitly temperature (or density)-dependent. On 
general thermodynamic grounds (see e.g. [35||) one can 
show that a is usually the first of the coefficients in the 
expansion @ that changes sign as temperature is low- 
ered. The sign change is a result of competition between 
either energy and entropy or different forms of entropy. 
The parameter T* accounts quantitatively for this com- 
petition and represents the spinodal temperature for the 
first order phase transition from the isotropic phase to 
the uniaxial nematic phase. As for the remaining param- 
eters: a > by definition and stability of the expansion 
requires e > and / > 0. Except for cases of multicrit- 
ical behavior, the signs of b, c, d, e, / are assumed not to 
change in the vicinity of T*. Hence, these coefficients, 
being weakly temperature-dependent, are assumed con- 
stants and taken at T = T*. 

The paramctrization of F in terms of q and u), Eq. (J5]), 
leads to a simple determination of absolute minima of 
F and, hence, a construction of the corresponding phase 
diagrams. Clearly, the form of Q, Eq. implies that 
the Iso-Nb and the Njj — Nb phase transitions can be 
either first- or second order. In other words we may ex- 
pect first-order, second-order and tricritical behavior at 
Iso-Nb and Nu — Nb transitions, depending on model 
parameters. 



III. PHASE DIAGRAMS 

Out of the five material parameters 6, c, d, e, / (cf> = 
f — e), introduced in Eq. (O, two are redundant and can 



be set equal to or ±1. This is a direct consequence 
of the freedom to choose a scale for the free energy and 
for Q. If not specified otherwise we choose e = 1 and 
c = 0,±1, and investigate the phase diagrams in the 
(a, 6)-plane as function of d and /. Additionally, we 
assume / > to guarantee the stability of the expansion 
(|6"|) against unlimited growth of q and replace / — e by 
4> (0 = / — e = f ~ 1) whenever convenient. We also 
make use of the free energy invariance with respect to the 
transformation: {6, d, q} — > {— 6, —d, —q}, which limits d 
to d > 0. The diagrams for d < are obtained as mirror 
images with respect to the 6 = line of those for d > 0, 
followed by a subsequent change of Nu± into Nu T . 

Interestingly, the relatively simple expansion © gener- 
ates a rich spectrum of possibilities for phase diagrams. 
We show that all of them can be divided into ten dis- 
tinct classes, where four involve only uniaxial phases. 
The remaining cases, corresponding to d < 0, are ob- 
tained from the classes discussed by applying the afore- 
mentioned 6 = mirror transformation. 



A. Phase diagrams with uniaxial phases: 

q ^ 0,u = 

We start by considering regions of stability of the uni- 
axial nematic. The necessary conditions for this phase to 
become, at least, locally stable read 



dF u 
dq 

d 2 F u 

dq 2 



q(a~bq + cq 2 + dq 3 + fq 4 







(9) 



= a-2bq + 'icq 2 + Adq 3 + 5fq 4 > 0. (10) 



The limit of local stability is attained when the inequal- 
ity (flTJ|) becomes equality, which, together with ([9]), de- 
scribes a saddle bifurcation in the model and represents 
spinodal lines. These conditions are particularly simple 
to solve for a and 6 in a parametric, q-dependent form. 
The nontrivial solution is 



cq 



2 1 2dq 3 + 3fq 4 



b = 2cq + 3dq 2 +4/g 3 , 



(11) 
(12) 



which, together with the trivial one: {a = 0, q = 0, (V6)} 
defines the borders of the area in the (a, 6) plane, where 
the solutions to the Eq. © are, at least, locally stable. 
Clearly, q runs over all real numbers. The subsequent 
calculation of the free energy at these local minima al- 
lows us to select the global minimum within the family 
of uniaxial solutions. 

A complete analysis of the model, including calcula- 
tion of the free energy, proceeds in a similar way. In par- 
ticular, we determine parametrically the transition line 
between the isotropic- and the uniaxial phases by solv- 
ing the system of equations: {F u = 0,dF u /dq = 0} for 
a(q) and 6(g). The solution reads 
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FIG. 2: Generic phase diagram with a direct, first order phase 
transition from isotropic- to uniaxial nematic phases, where 
gray color refers to phase transitions involving the Nu- phase. 
For b = a phase transition between isotropic and highly de- 
generated phase of q ^ but of arbitrary u) takes place at the 
quadruple point (T). For c > T became a tetracritical point. 
Parameters taken are (c, d, /) = ( — 1, 0, 1). Thin dashed lines, 
representing spinodal, are the solutions of the Eq. also 
a = spinodal is shown. 
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Subsequent analysis of the Eqs. (fTTI fl~3|) allows us to sin- 
gle out four topologically distinct classes of the phase 
diagrams with uniaxial- and isotropic phases. The rep- 
resentatives of each class are shown in Figs. [UGH The 
corresponding global stability sectors in the {c, d, /} pa- 
rameter space are given in Figs fTi]fT2l 

We shall now characterize each of the 'uniaxial' classes 
of the diagrams. 



1. Class (a) 

The first class is obtained for d = 0, Fig. [2] It contains 
a line of first order Iso *-* Nu+ phase transitions for 
b > 0, a line of first order Iso <-> Ny- phase transitions 
for b < 0, and a degenerated biaxial phase of q ^ and 
arbitrary ui, stable only along the 6 = line. We shall 
come back to the degenerated case in the last subsection 
of this paper. The Iso <-» Njj± lines have a common 
tangent a = 0. For c > the four lines meet at an 
isolated, tetracritical point, also often referred to as the 
Landau point. Its coordinates are (a, b) — (0,0). For 
c < the Landau point becomes a quadruple point of 
coordinates (a, b) = (||j,0), marked as 'T' in Fig. [5] 



FIG. 3: Generic phase diagram with a direct, Nu+ — Nu- 
first order phase transition. Parameters taken are c = d = 
/ = 1. For detailed definition of all symbols and lines used 
see caption to Fig. [2] The triple point (T) is localized at 
(a,b) = (0,0). 



The whole phase diagram is given analytically by 
_ y/((? + 16aff -c(c 2 - 48a/) 



16/ 



(15) 



where a > ^jG(c) with O being the step function. 

Now we concentrate on more complex cases with d ^= 
0. They are gathered in three classes of the diagrams, 
denoted (b)-(d). 



2. Class (b) 

Diagrams that belong to this class are similar to (a) 
except for an additional first-order phase transition line 
between the Njj+ and Nu- phases. A typical situation 
is shown in Fig. [3] Also in this case the transition line 
can be given in an analytical form as: 



a = u[2du(d + 2fu) + c{2d+5fu)]x 
[5c+2u{Ad + 5fu)] 
4(d + 5fu) 2 
u [3cd + 2u (3d 2 + 7fud + 5f 2 u 



2(d + 5/u) 



with the free energy 



F„ 



u 2 (3d + 5fu) [5c + 2u{4d + 5fu)\'' 
240(d+5/w) 3 



(16) 
(17) 

(18) 
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The parameter u must satisfy the inequalities 

<u< -j- for 25c/ <6c? 2 
5/ 5/ 



d 

-_<„<» 



for 25c/ > 6c! 2 



(19) 
(20) 



Generally, this topology is observed for the |£ parameter 
taken from outside of the interval {see discussion 



below leading to inequalities ([241 125)1 ) and is the most 
typical for the uniaxial family of the phase diagrams, 
Figs. [T21IT41 The appearance of the Nu+ — Nu- line is 
a result of competition between the third- and the fifth 
order invariants in the free energy expansion when the 
coefficients weighting these terms are of the opposite sign. 
For 25c/ > 9cP the three phases: Iso, N u+ and N V - 
meet at the triple point, T, of coordinates (a, 6) = (0, 0) 
where u = and F u — 0. At T the lines Iso — Njj + — and 
Iso — Njj- have a common tangent given by the a = 
line. For 25c/ < 6c! 2 the triple point moves away from 
the origin to a new location at 



3 (6c! 2 -25c/) 2 
a ~ 10000/ 3 
, 9c! (6c! 2 - 25c/) 
1 ~ 500/2 

which is obtained by substituting u 



(21) 
(22) 




FIG. 4: Generic phase diagram with a direct, Nu- — Njj- 
first order phase transition terminating at a critical point. 
Parameters taken are (c,d,f) — (1,1.05,0.306). Note that 
one of the spinodal lines connected with the Nu+ phase is 
separated from the transition line Iso — Nu+ by less than the 
thickness of the graph, and, hence, is invisible. For detailed 
definition of all symbols and lines used see caption to Fig. [2] 



-f^ into Eqs. (US]). 



3. Class (c) 



This class of the phase diagrams is perhaps the most 
interesting one among the uniaxial topologies. In addi- 
tion to the Njj + — Njj- transition line, shown in Fig. [3] 
and given by (|16|20j) . it also displays a direct Nu- — Nu- 
first-order phase transition line terminating at a critical 
point of the liquid-vapor type. 

Again this behavior results from the aforementioned 
competition between the third- and the fifth order terms 
in the free energy expansion. An example of the Nu- <-> 
Nu- line, together with the lines: Iso <-> Nu-, Iso <-> 
Nu+ and Iso *-* Nu+ <-> Nu- , is shown in Fig. [U The 
lines terminate at the Iso— Nu+ — Nu- triple point of co- 
ordinates (a, b) = (0, 0) and at the Iso—Nu- — Nu- triple 
point given by the formula (|2ip . The necessary condition 
for this class of the diagrams to appear is a requirement 
that the spinodal has two cuspidal points for the oblate 
states (q < 0). After inspecting q-dependence of the 
curve (fTTj) one easily finds that the (a, 6)-coordinates of 
these points are obtained by substituting 



Q± 



d±^d 2 -lcf 
4? 



(23) 



into Eq. (fTTj) . Additionally, the conditions < 8c/ < 
3c? 2 A c > (in our re-scaling c = 1) must be met 



for the cuspidal points with negative values of q to oc- 
cur. Taken together, these conditions guarantee that 
there exist two local minima (usually one of them be- 
comes the global one) and two local maxima in the free 
energy branch for the oblate states (q < 0). The local 
minima can finally convert into a stable Nu- — Nu- line, 
Eq. (HU), if 



6c! 2 < 25c/ < 9c! 2 



A 



c> 



d_ 

V 



J9d 2 - 24c/ 3c! 

— < u < . 

6/ " " 5/ 



(24) 
(25) 



Note that the conditions (|24I25|) are more restrictive than 
the ones for the cuspidal points of the spinodal. The 
first one, (|24|) . states that the Iso — Nu- — Nu- triple 
point disappears (and hence also the Nu- — Nu- line) 
for b > 0. Additionally, it guarantees the appearance of 
the Iso — Nu- — Nu- triple point (cusp) in the F u (q < 
0) = branch of the free energy for b < 0. The second 
inequality, (|2"B"|) , represents actually the same restrictions, 
but expressed in terms of u. Finally, coordinates of the 
critical point are obtained by substituting q-, Eq. (|23[) . 
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into Eq. l[TT]). This leads to 



>3/2 



4c/ (3d 2 - 2c/) - 3d 4 d (3d 2 - 8c/) ^ 



96/ 3 96V3/ 3 
9d 3 - 36c/d + y/3 (3d 2 - 8c/) 3/2 
W 2 • 



(27) 



Sector of stability of the class (c) is shown in Fig.rjJ] It is 
restricted to the area given by 6eP/25 < / < 9d 2 /25A/ < 
1. The richest phase sequence obtained for this class as 
temperature is lowered is I so — Njj + — Njj- — Njj—. 



4- Class (d) 

Quite interesting and untypical situation is met when 
cf/d 2 approaches one of its two limiting values in (f2"4"| . 
For cf = ^rd 2 , Fig. QH the Nu- - Njj- transition line 
and a part of the Njj + — Njj- transition line become 
reduced to a common straight line 



h= -2d a ' 



a e - 



8d 4 
625? 



■ 01 



(28) 



That is, the Nu- — Nu- transition line becomes also a 
line of triple points with the critical and triple point col- 
lapsing at a = — gffp! This case, illustrated in Fig. [5l 
makes us to expect that when higher orders in the ex- 
pansion ([S]) are taken into account the degeneracy of the 
Njj- — Njj- — Njj+ line should be removed and replaced 
by an (Iso) — Nu- — Nu- — Ntj+ bubble-shaped dia- 
gram with up-to three triple points. Bracket indicates 
that the branch Iso — Njj- does not need to be present. 
At cf = -^d 2 , which is the second of the two limits, the 
Njj- — Njj- line becomes reduced to a single critical point 
located at (a, b) = (gffp, 
the Iso — Njj- transition line. 

The phase diagrams described so far are stable against 
formation of the biaxial phase given that (/>=/ — 1 < 0. 
If this condition is fulfilled we can always find a uniaxial 
state with free energy lower than- or equal to the free 
energy of any biaxial state. Indeed, consider a biaxial 
phase of < u> < 1. A sufficient condition for the equi- 
librium value to b of lo in the biaxial phase is then given 
by 



wffi )■ The point belongs to 



OF 

duj 



F b (q) + -^q 6 LU b = 0, 0<co b <l, (29) 



which, when solved for F b (q) and substituted back to the 
biaxial free energy formula, Eq. (J5]), yields 

F Uax = F{q 7 u; b ) = F + F u (q) - \<f> q e W . (30) 

6 

The Eq. ([3D]) clearly shows that only for > (/ > 1) 
there is a chance to get a stable biaxial nematic phase. 
For <f> < the uniaxial state is always more favorable. 
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FIG. 5: Degenerated version of phase diagram shown in Fig. [4] 
The Nu- — Nu- line (together with the critical point) be- 
longs to the Nu- — Nu+ line starting at a = 0, both be- 
ing straight lines. The dashed straight line above the critical 
point is a continuation of the Nu- — Nu- line and serves 
as a reference to Nu- — Nu+ line. Parameters taken are 
(c,d,f) — (1,1,6/25). For detailed definition of all symbols 
and lines used see caption to Fig. [21 



The same conclusions are drawn for the biaxial state of 
u> = 1. By a direct calculation of the free energy we 
find for this case that the biaxial state of the free energy 
F(l = <lb,u = 1), is always less stable than one of the 
two uniaxial states: {q = ±q bl u> = 0}, where q b is value 
of q in the biaxial phase. 



B. Phase diagrams with biaxial nematic phase: 

<7^0,0<u<l 

The discussion of the previous section shows that, 
generally, a stable biaxial nematic phase is found for 
4> = f — e = f — 1>0. In this section we analyze this 
case more thoroughly. We start by pointing out that the 
sign of the F b (q) term in Eq. ^ decides about the rela- 
tive stability of the biaxial order with respect to all other 
phases involved. Generally, a uniaxial phase becomes un- 
stable against formation of the long-range biaxial order 
if F b (q) < 0, which implies that 



b^dq 2 



<f>q d 



q^Q. 



In addition, for F b (q) < —^(f>q e , or, equivalently 



b^-dq 2 



q^0 



(31) 



(32) 



the phase biaxiality, u>, attains its maximal value u> b = 1, 
Eq. (|29[) . The equality sign in the condition ([3"Tj) marks 
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a bifurcation from the uniaxial to the biaxial phase. To- 
gether with © , this can be solved for a and b to give the 
spinodal lines in a parametric form: 



(a,b) 



-cq 2 -^q 3 
5 



4 3d 2 
5 



(33) 



A few general conclusions can be drawn from the for- 
mulas (pl3^1> and inequality (jnH- First of all, if Eq. ^ 
is fulfilled on a globally stable uniaxial nematic branch, 
the transition Njj — Nb is second order. Satisfying rela- 
tion (|33p on a locally stable uniaxial branch results in a 
first order Njj (I so) — Nb phase transition. That is, the 
bifurcation scenario allows for a possibility of a tricritical 
point on the Nu — Nb line. Second order I so — Nb tran- 
sition is only admitted to states of maximal biaxiality 
( W = l). 

For the biaxial branch of the free energy a more quan- 
titative analysis can be given. In particular, the biaxial 
free energy (|30|) can be expressed in an equivalent form 
as 



F, 



biax 



&2 

60 



1 



1 



-aq 



7lQ 



(34) 



55b 



25a 2 (q 2 + c) 0] (38) 
6d 2 g 2 ™250 [q*-ql (d + f qi ) + 

c(q 2 -q 2 )] (39) 
5O<? 2 = Ad 2 - 400ft d - 5OC0 ± 

\f2~s/5 2 {2d 2 + 4O0«id + 25(2C - c)0), (40) 

with S = 2d + 50ai and £ = jq\ + c. The parameter qi 
runs over the uniaxial branches, where q\ > for Njj+ 
and qi < for Njj_. Analyzing various cases we are 
able to single out six additional classes of the diagrams, 
shown in Figs. IBI fTTl that supplement the uniaxial family. 
Again, the phase diagrams with the Nb phase should be 
correlated with Figs. fT2lfT4l where sectors of absolute sta- 
bility of a given class are shown in the {c, d, /} parameter 
space. 

The discussion of the biaxial phase diagrams will pro- 
ceed in a similar way as for the uniaxial case, that is 
we again start with the case of d = 0. In this limit we 
can distinguish between the two different classes of the 
diagrams, all being symmetric with respect to the b = 
line. 



where 



2bd 

50 ' 



7 



6d 2 
250 



and 7 2 > 4a. (35) 



A convenient parametric form for the I so — Nb line now 
easily follows from the equation F^ax = 0, supplemented 
with the condition for qi,: (dFt > i ax /dq) q=qb = 0. The 
solution of practical importance may be expressed as 



2bd 

' 50 



79 



o\b 2 



3 7 )0 



(36) 

Additionally, a stability criterion of the biaxial solution is 
given by the condition that determinant of second deriva- 
tives of the free energy is positive definite. This means 
that the biaxial phase is locally stable if 



d 2 Fd 2 F 



d 2 F 



> 



dq 2 dui 2 \dqdu> 

4a + 7 ^^/t 2 ~ 4a — 7^ 



> 0. 



(37) 



The limiting case of vanishing determinant gives two 
straight lines in the (a, 6)-plane: {a = 0, 7 2 = 4a}, which 
are further spinodals of the model. 

Detailed analysis of relative stability of Iso, Njj and 
Nb phases shows that all 'uniaxial' phase diagrams, 
Figs. [HO have their biaxial counterparts. Generally, the 
biaxial phase replaces, at least partly, the Njj+ — Njj- 
transition line by the two lines: Njj + <-> Nb and 
Nu- Nb- They can be given in a parametric form 
as functions of the real parameter q\ 

55a = qi [(6q 2 - lOqj) d 2 - 10( qi d+ 



1. Class (e) 

The first class, shown in Fig. [SI is similar to (a), Fig. [5] 
The only difference is that the line separating Njj+ and 
A^- splits itself into Njj + <-> Nb and Nu- Nb lines 
of second order phase transitions with JV^ phase posi- 
tioned in between. We find this class stable for c > 
and / > 1(= > 0). As previously the uniaxial lines are 
given by (fT"5|) . For the Njj — Nb lines the formulas 
now simplify to 



b 2 



i (3ac - c 3 + Vc 2 - 4a (c 2 - a)) 2 



(41) 



where a < 0. The four phases: Iso, Njj + , Nb and Njj- 
meet at the Landau (tetracritical) point: L = (a, b) = 
(0, 0). Additionally, for c = the N u+ - N B and N v _ - 
Nb lines have a common tangent at L, which is given by 
the a = 0-line. For c = 1 this tangent is the b = 0-line. 



2. Class (f) 

The diagrams of this class, Fig. [3 are also derived from 
(a) and observed when c < 0. Again the uniaxial lines 
are given by Eq. (15]) and the N v - Nb lines by Eq. ([?!) . 
The latter represent thermodynamically stable, second- 
order transition lines if 



a < 



a < 



3c 2 (2/ - 1) 
4(/ + l) 2 



for 1 < / < 2 (42) 
for / > 2. (43) 
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FIG. 6: Generic phase diagram for d = and c > 0. Biaxial 
phase is sandwiched between two uniaxial phases. A direct 
phase transition from the isotropic phase to biaxial nematic 
is possible through the Landau point. Solid lines represent 
phase transitions of first order, dashed lines second order. 
Two cases are shown: c = 0,d — 0, / = 1.5 (thin lines) and 
c = l,d — 0, / = 1.5 (thick lines). As previously, gray lines 
represent phase transitions involving Nu- phase. 



A new feature shown is a splitting of the Landau point 
into two triple points where Iso, Njj and Nb meet. The 
position of the triple points is 



(a,b 2 ) 



3c 2 (2/-l) 27c 3 (/ - If 
4(/ + l) 2 ' 8(/ + l)3 



(44) 



Both triple points are connected by a direct Iso «-> Nb 
line of first order phase transitions for which we have 



b 2 = -6 



(c 2 - 6a) - (c 2 - 4a) 



3/2 



(45) 



Depending on /, the phase transition between Njj and 
Nb can be either first or second order. For 1 < / < 2 
only second order Njj — Nb transitions are realized. For 
/ > 2 the second order transition line (|41[) is separated 
from the triple point by the Njj — Nb line of first order 
phase transitions (u / 0). Both Njj — Nb lines meet at 
the tricritical point given by 



2\ _ I c „3i2 



( a y) = [ -,--c 



(46) 



FIG. 7: Generic phase diagram for d = and c < 0. Biaxial 
phase becomes stable between two uniaxial phases. A di- 
rect phase transition from isotropic phase to biaxial nematic 
is possible along the line between the two triple points, T, 
that replace the Landau point. Solid lines represent phase 
transitions of first order, dashed lines second order. Tcp 
stands for tricritical point. Two cases with / «s 2 are shown: 
c = -1, d = 0, / = 1.7 (thin lines) and c = 1, d = 0, / = 3.1 
(thick lines). For the meaning of lines see caption to Fig. [6] 



the diagrams that emerge in this case by referring directly 
to the parameter space division as shown in Figs. [L21Tl"4l 
New classes will be parameterized by c = ±1, 0. It turns 
out that for c = — 1, the effect of nonzero d is merely to 
distort the phase diagrams classified as (f). The distorted 
diagrams that preserve all features of (f) are separated 
from the new class (f), Fig. [51 by curves: 



rf2 _ 25c (v^-1) 



6(^.-2^-2) 



< <f> < 4 + 2^3. (47) 



The curves are pictured dark-gray in Fig. 1121 The area 
to the right, shown as light-gray, represents the class (f). 
The class differs from the deformed versions of (f-Tcp)- 
like diagrams with two tricritical points and of (f) with- 
out tricritical points by the presence of one tricritical 
point on the Nu + — Nb transition line. 



4- Class (g) 



3. Class (f) 

Now we turn to a more complex case, namely that of 
d ^ 0. It is quite convenient to discuss new features of 



The case c = 0, Fig. [T31 results in a new class of the 
diagrams shown in Fig. [9] One of the differences between 
(g) and (P), exemplified in Fig. [51 is the absence of the 
direct transition between Iso and Nb- The Nb phase 
branches off the Njj + — Njj- first-order transition line 
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FIG. 8: Generic phase diagram for d / and c < 0. Biaxial 
phase becomes stable between two uniaxial phases. A direct 
phase transition from the isotropic phase to biaxial nematic 
is possible along the line between two triple points (T) that 
replace the Landau point. Only one tricritical point (Tcp), 
along Nu+ — Nb, is possible for this class of diagrams. Pa- 
rameters taken are: c = —1, d = 1, / = 3. For the meaning of 
lines see caption to Fig. [5] 



at the Njj+ — Nb — Nu- triple point. Interestingly, we 
observe a maximum along the second-order Nb — Njj- 
transition line at the location given by 



(a, b) = 



16^(6 - 1) 4d 3 



6250 4 



125</> 2 



(48) 



This maximum indicates that we can observe reentrant 
biaxial nematic phase as temperature is lowered. Conse- 
quently, it leads to a very rich sequence of phase transi- 
tions, e.g. I so — Nu+ — Nb — Njj^ — Nb- The reentrant 
phase and hence also the maximum disappear for / > 2. 
In the interval 2 < / < 1 + 2/V3, shown as sector (g') 
in Fig. 1131 the remaining features of the diagram, Fig. [51 
are left unchanged. 



5. Class (h) 

For c = 1 we identify two new classes of the diagrams, 
denoted (h) and (i). The class (h), Fig. [lOj is derived 
from (e), the difference again being the presence of max- 
imum along the second-order Nb — Njj- transition line at 
(a, b) given by Eq. (|48| . That is we again can observe a 
sequence of phases with reentrant biaxial nematic. Sector 
(h), Fig. |T?1 is separated from the neighboring sectors (g) 
and (h+c) by the following lines: the dashed one given 



FIG. 9: Generic phase diagram for d 7^ and c = 0. Biaxial 
phase becomes stable between two uniaxial phases. A direct 
phase transition from isotropic phase to biaxial nematic is not 
possible. Two triple points (T) are connected by the Nu+ — 
Nu~ line of first order phase transitions. One tricritical point 
(Tcp) appears on the Njj+~Nb line. Maximum along Nu- — 
Nb allows for two biaxial nematic phases on the temperature 
scale, separated by the Nu- phase. The low-temperature 
biaxial phase is referred to as reentrant Nb ■ Parameters taken 
are: c = 0, d = 2, / = 1.6. For the meaning of lines see caption 
to Fig. H 



by / = 1 + 6d 2 /25 (0 < d < 5/VS) and the continuous 
one given by / = 9d 2 /25 (d > 5/^3). 



6. Class (i) 

This class of the diagrams is essentially a combination 
of (g) and (c) and yields the richest sequences of phases 
and of corresponding phase transitions. They include 
reentrant biaxial nematic, Njj+ — Nb tricritical point 
and a line of phase transitions between identical uniaxial 
phases terminating at a critical point. Exemplary phase 
diagram is given in Fig. [TTJ Sector of stability for this 
class, denoted (i) in Fig. [TH is limited by the follow- 
ing curves: / = 1 + 6d 2 /25 (d > 5/V3), / = 9d 2 /25 
(5/3 < d < 5/V3), / = 1 (5/3 < d < 5/V&) and 
/ = 6d 2 /25 (d > 5/v6)- In a small sector, named (h+c), 
the tricritical point and the Nu + — Njj- line disappear, 
the resulting phase diagrams being a combination of (h) 
and (c). 
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a 

FIG. 10: Generic phase diagram for d / and c = 1. This 
class of the diagrams is a deformed versions of (e). A major 
difference between (e) and (h) is a maximum along Nu- — Nb 
that allows for reentrant Nb- Parameters taken are: c = 
1, d = l,f = 1.5. For the meaning of lines and of reentrant 
Nb see captions to Figs. [6] and [9] respectively. 




^(a) (b) 

1 2 
d 



FIG. 12: Sectors, in {c = —1, d, /} parameter space, of abso- 
lute stability of phase diagrams labeled from (a) through (i). 
(f-Tcp) stands for (f)-class with two tricritical points, Fig. [7] 
Deformed versions of (f) and (f-Tcp) diagrams are realized 
within sectors marked white. 
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FIG. 11: Generic phase diagram for d / and c = 1. It 
combines properties of (g) and (c). Parameters taken are: 
c = l,d = 3,/ = 2.75. 



FIG. 13: Sectors, in {c = 0, d, /} parameter space, of abso- 
lute stability of phase diagrams labeled from (a) through (i). 
(g') stands for (g)-class without reentrant Nb (i.e. without a 
maximum along Nu- — Nb)- 
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d 

FIG. 14: Sectors, in {c = l,d, /} parameter space, of abso- 
lute stability of phase diagrams labeled from (a) through (i). 
(h+c) stands for diagrams that combine properties of (h) and 
(c) (see Figs. HTTOjl . 



C. Degenerated phase diagrams for b=0 

This case requires a few comments for not all phase 
sequences with b=0 can be identified from the diagrams 
that we have given so far. There are also subtle symmetry 
issues due to the presence of accidental degeneracy. 

The identification of phases and of phase sequences for 
b = is straightforward for d =/= 0, in which case they 
follow directly from the diagrams representing classes 
(b,c,d,f',g,h). Only the case d = 0, represented by the 
diagrams (a,e,f), requires separate comments. By in- 
specting the free energy expansion (J6]) for b = d = we 
find immediately that only three states can be realized at 
equilibrium: (i) degenerated uniaxial state of u> = for 
/ < 1; (ii) degenerated biaxial-uniaxial state of arbitrary 
lu for / = 1, and (Hi) biaxial state of maximal biaxiality 
(u> = 1) for / > 0. 

The uniaxial state, denoted (i), has the same free en- 
ergy for oblate and prolate states, which means that the 
system creates oblate and prolate domains with the same 
energy cost. Due to this accidental degeneracy the sym- 
metry group of the state can be classified as 2?ooh x Z-x . 
Consequently, the transition from the isotropic phase to 
the degenerated uniaxial phase can be either second- 
order (c = 1) or first order (c = —1) with a tricritical 
point at c = 0. The transition temperature and the or- 
der parameter for temperatures below the transition are 
given by {a — for c = 0, 1; a = y|j for c — — 1} and 

q 2 = ^ C ~~^^ — , respectively. 



The states (ii,iii) have mathematically the same form 
of the free energy as one for the degenerated uniaxial 
case. All formulas are reproduced from the case (i) if we 
substitute / = 1 there. Hence, again, the phase tran- 
sition from the isotropic phase to the corresponding or- 
dered phase can be either first or second order with an 
intermediate tricritical point. The difference between the 
cases (ii) and (Hi) is in symmetry. For (ii) the only re- 
striction on Q is Tr(Q ) = const. The accidental degen- 
eracy of the equilibrium solutions for Q is the rotational 
invariance in five-dimensional space of components Q a p. 
The relevant symmetry group is thus 0(5), which seizes 
up the difference between uniaxial and biaxial domains. 
We call it the degenerated biaxial-uniaxial phase. 

In the case (Hi) the tensor Q is given by Eq. (P) with 
e.g. qo = 0, that is by one of the three degenerated solu- 
tions fulfilling maximal biaxiality condition Tr(Q 3 ) = 0. 



IV. DISCUSSION 

The Landau-deGennes theory of biaxial nematics pre- 
sented in this paper has been elaborated to show the out- 
come of mathematical structure of the expansion that is 
based solely on the order parameter Q. The phenomcno- 
logical approach is particularly simple and the knowledge 
of full spectrum of predictions of one of the most com- 
monly cited theories is desirable, particularly because of 
current interest in seeking for stable thcrmotropic biax- 
ial nematic phase. The analytical formulas are given for 
almost all transition lines, characteristic points of the 
lines, and for the stability range of a given class of the 
phase diagrams. Except for the purely uniaxial group of 
the diagrams for / < 1, the biaxial phase is naturally 
stabilized between prolate and oblate uniaxial nematics. 
Phase transitions to the biaxial phase can be cither first 
or second order with a possibility of a tricritical point. 
Due to a competition between cubic and fifth-order in- 
variants the direct Nu_ — Njj-, Nu + — Njj- and the reen- 
trant biaxial nematic phase are also possible. The Lan- 
dau I so — Njj+ — Nb — Nu- (tetracritical) point can split 
into two triple points positioned either on the I so — Nb 
transition line or on the Njj+ — Nu_ line. 

One may wonder then why the biaxial nematic phase is 
so difficult to account for experimentally. The practical 
difficulty could be that for real systems the parameter 
/, responsible for the stabilization of the biaxial nematic 
phase, is much too small compared to other coefficients of 
the Landau expansion, so we effectively stay in the uniax- 
ial sector of the diagrams (/ < 1). An alternative expla- 
nation could be that smectic and crystalline phases, not 
taken into account, may interfere before the right ther- 
modynamic parameters are reached. Clearly, the best 
choice of the Landau coefficients to get Nb absolutely 
stable would be that where Nb bifurcates directly from 
the isotropic phase. At microscopic level it would then be 
of interest to construct molecular models showing generic 
types of the diagrams identified phcnomcnologically. To 
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this goal it is necessary to establish a bridge between 
molecular and phenomenological approaches, in particu- 
lar one needs a molecular interpretation of the alignment 
tensor and, at least, of the (a, 6)-parameters entering 
the expansion. The problem is relatively simple in the 
mean-field approximation and the solution has already 
been given [lfl [3f| for the class of the so called L = 2 
models with 2?27i-symmetric hard molecules/soft interac- 
tions [H,[2l||. A PPly m g formulas (14-20) from J^lto the 
mean-field versions of the models [l?], [20 [2ll |22| l36l HtJ 
we recover diagrams represented by (e) [13, H2, \2M > (f ) 
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